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Abstract 

A quantum mechanics analogy is used to determine the forces act- 
ing on and the energies of solitons governed by the nonlinear Schrodin- 
ger equation in finite intervals with periodic and with homogeneous 
Dirichlet, Neumann and Robin boundary conditions. It is shown that 
the energy densities remain nearly constant for periodic, while they 
undergo large variations for homogeneous, boundary conditions. The 
largest variations in the force and energy densities occur for the Neu- 
mann boundary conditions, but, for all the boundary conditions con- 
sidered in this paper, the magnitudes of these forces and energies 
recover their values prior to the interaction of the soliton with the 
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boundary, after the soliton rebound process is completed. It is also 
shown that the quantum momentum changes sign but recovers its orig- 
inal value after the collision of the soliton with the boundaries. The 
asymmetry of the Robin boundary conditions shows different dynamic 
behaviour at the left and right boundaries of the finite interval. 

KEYWORDS: Nonlinear Schrodinger equation, two-point initial- 
value problems, quantum analogy, forces 

1 Introduction 

The one-dimensional nonlinear Schrodinger equation is ubiquitous in many 
branches of mathematics, physics and engineering [T] and [2], e.g., fluid 
mechanics, plasma physics, nonlinear optics, nonlinear wave propagation, etc. 
Its initial- value problem is integrable by the inverse scattering transform and 
has a truly remarkable analytical structure. Despite the great interest that 
this equation has received in the past two decades through analytical and 
numerical studies, the existence and uniqueness of the nonlinear Schrodinger 
equation in the quarterplane, i.e., for semi-infinite spatial domains, seems to 
remain an issue of current research [3]. 

Three main lines of research on the nonlinear Schrodinger equation in the 
quarterplane have been undertaken. The first is based on results obtained 
numerically and uses ad hoc assumptions [3J. The second one is based on the 
use of the inverse scattering transform and Backlund transformations [5], 
while the third one employs the infinite sine Fourier transform [6] . However, 
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these three lines of research have been concerned with the generation of 
solitons at the boundaries, e.g., boundary-generated solitons. 

Ablowitz and Segur [7] showed that the quarter-plane, initial-boundary 
value problem is solvable by means of the nonlinear, sine Fourier transform 
if certain symmetry conditions hold. In particular, Ablowitz and Segur [7j 
proved the solvability of semi-infinite line problems subject to homogeneous 
Dirichlet or Neumann conditions. Fokas [H] and Bikbaev and Tarasov [5] have 
extended this result to semi- infinite line problems subject to homogeneous 
Robin boundary conditions. Fokas [6] has studied the generation of solitons 
by a non-homogeneous Robin boundary condition in the semi-infinite line 
with zero initial conditions. 

Some authors believe in the almost integrability of the nonlinear Schro- 
dinger equation in the quarterplane, also called 'forced integrable nonlinear 
Schrodinger equation' [U] because the boundary conditions can be viewed as 
forces that act at the boundaries. Kaup [5] showed that the difficulty in the 
solvability of the forced Dirichlet initial-boundary value problem in the quar- 
terplane lies in the determination of the spatial derivative of the amplitude 
at the boundary, and in the determination of its primitive in the Neumann 
case. However, his studies have not been able to provide a connection, if 
there is any, between solitons in initial-value and quarterplane problems. 
Other authors, such as Kaup [8] and Kaup and Hansen jl] have studied 
both analytically and numerically the generation of solitons in quarter-plane 
problems subject to non-homogeneous Dirichlet boundary conditions. 

If little is known about the existence and uniqueness of solutions of the 
nonlinear Schrodinger equation in the quarterplane, much less is known about 
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that equation in finite intervals, i.e., two-point, initial-boundary value prob- 
lems, because neither the inverse scattering transform nor the infinite sine 
Fourier transform can be applied to them. Numerical solutions, however, may 
be obtained for these problems by employing, for example, global spectral, 
finite difference or finite element methods. In this paper, a numerical study 
of the nonlinear Schrodinger equation subject to periodic and to homoge- 
neous Dirichlet, Neumann and Robin boundary conditions in finite intervals 
is presented in order to determine the forces on and the energies of the soli- 
tons by using a quantum mechanics analogy. An analogy similar to the one 
presented in this paper has been previously used by Bian and Chan [10J who 
only studied the nonlinear and diffraction forces on solitons governed by the 
initial-value problem of the nonlinear Schrodinger equation. These authors, 
however, used a filter and their force and energy densities are not consistent 
with a strict quantum mechanics analogy. 

2 The One-Dimensional, Nonlinear Schrodin- 
ger Equation 

The one-dimensional, time-dependent, nonlinear Schrodinger equation with 
a cubic potential can be written as 
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where the potential well has the following form 



V = -q— \u\ 2 
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t is time, x is the spatial coordinate, u is a complex amplitude, % is the 
Planck constant, q is a coefficient, m is the mass of the soliton considered as 
a quantum particle, and the subscripts denote partial differentiation. 

A quantum-like analog of the nonlinear Schrodinger equation can be eas- 
ily obtained by assuming that the complex amplitude, it, of the nonlinear 
Schrodinger equation is like the wave function of a quantum particle. This 
allows to determine its quantum momentum, quantum energy and forces as 
follows. The linear quantum momentum is defined by the one-dimensional 
operator [11] 

< (3) 
whose mean value is given by the following time-dependent expression 

(P) = %M = - ih ff***f (4) 
\u,u) \ v u*uax 

where V is the spatial domain, and it* denotes the complex conjugate of it. 
The quantum momentum density may be defined as 

p(x,t) = —ihu*Ux- (5) 

The linear quantum energy is defined by the operator 

d 

i— 
dt 

whose mean value is 



d 

E = ih—. (6) 



(E) = %^ = ih Sf^f . (7) 
{u, u) j v u*u dx 

The energy density is 

e(x, t) = ihu*Uj:. (8) 
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Equation ([!]) can be nondimensionalized using the linear transformations 



pm - / h 

u — 7-J— w, t — > t, x — > \ — x (9) 
y qh V 2m 

to yield 

iut = —u xx — \u\ 2 u (10) 

Hereon, unless otherwise stated, nondimensional variables will be used. 
Note that the above dimensional definitions of the quantum momentum and 
energy may be used to obtain their dimensionless counterparts by simply 
setting h and m equal to unity in the definitions of the dimensional quantum 
momentum and energy densities. 

The energy density can be written as 

e(x, t) = -{u*u xx + H 4 ) = e k (x, t) + e v (x, t) (11) 

where the kinetic and potential energy densities are, respectively, 

e k (x,t) = -u*u xx , e v (x,t) = -|w| 4 . (12) 

The spatial derivatives of the energy densities define effective force den- 
sities as follows. From the potential energy density, the following nonlinear 
force density that produces self-focusing on solitons is obtained 

d d 

fn(x,t) = -~e v (x,t) = — (M 4 ) (13) 

while, from the kinetic one, the following diffraction force density that is 
responsible for the diffraction effect on the soliton results 

f d (x, t) = -—e k (x, t) = — (u*u xx ) . (14) 



The kinetic energy and diffraction force densities defined above have com- 
plex values. In quantum mechanics, only the mean values of the energy and 
force can be physically measured and these values are always real due to the 
Hermitian quantum operators. 

In order to obtain physical insight from the energy and force densities 
defined above, it is convenient to use real values by taking either their real or 
imaginary part. In this paper, the real parts of the energy and force densities 
are used. 

A physically and mathematically consistent use of the quantum analogy 
may also be obtained by defining the energy density as e/g where g is a 
constant. For example, the following renormalization factor may be used for 
the initial value problem of the nonlinear Schrodinger equation 

g — u*u dx (15) 
Jv 

which coincides with the first invariant of the initial-value problem of the 
nonlinear Schrodinger equation. However, this value of g may not be constant 
for two-point, initial-boundary- value problems (cf. Section 4), i.e., it may 
not be used to obtain a physically and mathematically consistent quantum 
analogy for two-point, initial value problems. 

Bian and Chan [10] defined the energy density as e/g where g=u*u in 
order to assess physically that a soliton can be considered as the result of 
two opposite phenomena, i.e., the diffraction and the self-focusing caused 
by dispersion and nonlinearity, respectively, for the initial-value problem of 



Eq. (10). Since their renormalization factor, g, is a function of both space 
and time, its use is not consistent with a quantum mechanics analogy for 
the nonlinear Schrodinger equation. As a consequence, the energy and force 



densities obtained by Bian and Chan are somewhat artificial, since their filter 
or renormalization factor affects the values of these densities in a different 
manner depending on the soliton location at each instant of time. 



The 1-soliton solution of the initial- value problem for Eq. (10) is 



u(x,t) = A sech(£) exp(w/). 



(16) 



where 



£ = A^{x-x -ct) 

1 c 2 
V = 2 \c(x - x ) + (A 2 - -)t + <j) ], 



(17) 
(18) 



A, c, xo and (fio are the soliton's amplitude, velocity, initial position and 
initial phase, respectively, and the nonlinear and diffraction force densities 
for g—1 are, respectively, 



/„ = -2V2A 5 sech 4 £tanh£ 



(19) 



A 3 

fd = -£ sech 2 £ [-QAci sech 2 £ + 8V2A 2 sech 2 £ tanh £ + 
+ V2(c 2 - 2A 2 ) tanh£ + AAci 



(20) 



while the total force density, i.e., the sum of the nonlinear and diffraction 
force densities, is 



A 3 

ft = ~£ sech 2 £ [-6Aci sech 2 £ + V2{c 2 - 2A 2 ) tanh£ + AAci 



(21) 



The above equations show that the nonlinear force density is less than 
zero behind the location of the soliton's maximum amplitude and greater 
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than zero in front of it, whereas the real parts of both the diffraction and 
total force densities depend on the values of c and A. 

The numerical results of Bian and Chan [10] indicate that their total 
force density, i.e., the sum of the real part of the diffraction force and the 



nonlinear force densities, vanishes for the 1-soliton solution, cf. Eq. (16). For 
the exact iV-soliton solution [131 E] > the energy density of Bian and Chan 
gives a nonvanishing force dominated by the nonlinear density force which 
indicates that the solitons are compressed while propagating. However, their 
interpretation can be criticized on the grounds that diffraction effects on the 
soliton must also be considered and that their renormalization factor which is 
a function of both space and time, affects in different manner the local force 
densities. Furthermore, it is important to note that their definition can cause 
computational problems when calculating the nonlinear and diffraction force 
densities due the small denominator introduced by the filter in their energy 
expression, especially in the soliton tails. 



3 Boundary Conditions 

In this paper, the nonlinear Schrodinger equation is studied numerically in a 
symmetric, finite interval T>=[—L,L] subject to the following homogeneous 
boundary conditions 

u(-L,t)+ju x (-L,t) = 0, u(L,t) + ju x (L,t) = and t > 0. (22) 

The values 7=0 and oo correspond to Dirichlet and Neumann boundary 
conditions, respectively. 
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Since the nonlinear Schrodinger equation and the homogeneous Dirichlet 
and Neumann boundary conditions are invariant under mirror reflections in 
x, their respective initial-boundary-value problems also exhibit this invari- 
ance; as a consequence, the interaction of a soliton with the left boundary 
is identical to that with the right boundary. This symmetry or invariance is 
lost if 7 is finite and different from zero, i.e., if homogeneous Robin boundary 
conditions apply at both boundaries. In this paper, mixed boundary condi- 
tions corresponding to 7=1 are considered, and the results for the Dirichlet, 
Neumann and Robin boundary conditions are compared with those for both 
the initial-value problem and the periodic boundary conditions 

— (x,t) = —(x+2kL,t), Vn>0, keZ, xeP=[-L,L], t > 0. 

(23) 

4 Relationships between the Quantum Mo- 
mentum and Energy Densities 

The most important invariants of the initial-value problem of the nonlinear 
Schrodinger equation are the first, second and third ones. The first, known 
as wave mass or 'number of particles', is 

/ \u\ 2 dx, (24) 
Jv 

whose integrand will be denoted by pn- The second invariant represents the 
total momentum in the Hamiltonian formalism and is given by 

/ % (u*u x — uu* x ) dx (25) 
Jv 
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2 --\u\ A )dx (26) 



and whose integrand will be referred to as p/ 2 . The third invariant is the 
total energy or Hamiltonian, i.e., 

1 

/dV - * 1 2' 

whose integrand is the Hamiltonian density of the nonlinear Schrodinger 
equation and will be denoted by p/ 3 . 

The real and imaginary parts of the quantum momentum density are 

i 

Ke{p(x, t)} =p R = -- (u*u x - uu* x ) (27) 

lm{p(x, t)} = pi = -- {u*u x + uu* x ) (28) 
those of the kinetic energy density 

ne{e k {x, t)} = e R = ~ {u*u xx + uu* xx ) (29) 

lm{e k (x, t)} = e z = - (u*u xx - uu* xx ) (30) 

and those of the diffraction force density 

1 d 

TZe{f d (x,t)} = f R = -— (u*u xx + uu* xx ) (31) 

i d 

lm{f d (x,t)} = // = --— (u*u xx - uu* xx ) (32) 

From the above equations, the following relations may be easily deduced 
d o 0p R dp I2 1 d 

m pn = - 2 -te = -te> PI = -2te pn - (33) 

It can also be easily shown using the nonlinear Schrodinger equation that 
the following equations hold 

dpR dp n , , dp 



dt dx ' dt 
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+ fR, = fi, (34) 



dp dp I3 de k dp I3 9 

m=-^-^ -W = d- x ^ Ut + u ^ ] (35) 

For initial-value problems, i.e., T>=(— 00,00) subject to \u\ — > as \x\ — > 

00, the above equations imply that the total mass or number of particles, 
the total momentum in the Hamiltonian formalism and the total energy are 
invariant. Furthermore, 

/ pdx = 0, (36) 

1. e., the total quantum momentum is zero. 

For quarter-plane or semi-infinite, initial-boundary- value problems, i.e., 
T>=[0, 00) subject to \u\ — > as x — > 00, the total number of particles 
is constant for homogeneous Dirichlet or homogeneous Neumann boundary 
conditions at x — 0, and the imaginary part of the total quantum momentum 
is only zero for homogeneous Dirichlet boundary conditions at x = 0, whereas 
both the total quantum momentum and its real part are different from zero. 

For the finite-line, initial-boundary-value problems considered in Section 
3, the total number of particles is constant for homogeneous Dirichlet or ho- 
mogeneous Neumann boundary conditions at both boundaries; the imaginary 
part of the total quantum momentum is zero for periodic or for homogeneous 
Dirichlet boundary conditions at both boundaries. The total quantum mo- 
mentum and its real part are different from zero. 

5 Numerical Method 



The Crank- Nicolson method was used to discretize Eq. (10) in the interval 
[— L, L], and a Newton- Raphson technique was employed to solve the re- 
sulting system of nonlinear algebraic equations. The diagonally dominant, 
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tridiagonal system of linear algebraic equations that result from the Newton- 
Raphson method was solved by means of a 2 x 2 block-oriented version of the 
Thomas algorithm. In the periodic natural optimization of the Gaus- 

sian elimination technique for periodic problems that yields quasi-tridiagonal 
systems has been used before applying the Thomas algorithm [12] . 

The Crank-Nicolson method is conservative since it preserves a discrete 



equivalent of the first invariant, cf. Eq. (24), i.e., the L norm of the solution, 
for the discrete, initial-value problem, and for the two-point initial-value 
problem with periodic or Dirichlet boundary conditions. 

The evaluation of the high-order derivatives of the soliton amplitude re- 
quired for the calculation of the soliton energy and force densities was per- 
formed so as to preserve both the maximum possible symmetry in the finite 
difference operators and second-order accuracy. It should be noted that, in 
the periodic case, symmetric, second-order accurate in space, stencils were 
employed at all the spatial grid points, while, for the Dirichlet, Neumann 
and Robin boundary conditions, asymmetric stencils were introduced near 
the boundary points in order to avoid the use of fictitious points. 



6 Presentation of Results 

In this section, some sample results (cf. Figs. 1-5) that illustrate the non- 
linear force density, the real part of the diffraction force density, the total 
force density, i.e., the sum of the nonlinear force density and the real part 
of the diffraction force density, and the real part of the momentum density 
are presented as functions of space and time for the four types of boundary 
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conditions considered in this paper. Figures 1-5 also show the space-time 
isocontours of the three-dimensional data presented in these figures. The re- 
sults presented in Figs. 1-5 as well as in Figs. 6-10 correspond to an interval 
of L=50, A—c—1, x o =0o=O, and spatial and temporal step sizes equal to 
0.25 and 0.01, respectively. 

Figure 1 corresponds to periodic boundary conditions and shows that the 
nonlinear force density has an S-shape which is nearly the mirror reflection 
of that for the diffraction force density. The total force density has also an S- 
shape similar to the nonlinear one which indicates that, for periodic boundary 
conditions, the nonlinear force density is larger than the diffraction one. The 
momentum density has a bell-shape similar to that of the soliton amplitude. 

Figure 2 illustrates the large changes in the magnitude of the force den- 
sities introduced by homogeneous Dirichlet boundary conditions and, espe- 
cially, in the collision of the soliton with the right boundary. It is interesting 
to note that the shape of the nonlinear force density prior to and after the 
collision of the soliton with the right boundary is the same owing to the 
symmetry of the nonlinear Schrodinger equation and boundary conditions; 
however, its magnitude increases greatly in the collision process because the 
amplitude of the soliton at the boundary is zero. Figure 2 also shows that 
the nonlinear force density is larger than the diffraction one, and that the 
momentum density changes sign once the soliton rebounds from the right 
boundary. Note that the Dirichlet boundary conditions require that the 
momentum density at the boundary be zero; therefore, the change in the 
momentum is smooth and the change in its sign is associated with the dif- 
ferent directions of the soliton velocity prior to the collision and after the 
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soliton rebounds from the boundary. 

Figure 3 corresponds to homogeneous Neumann boundary conditions and 
illustrates the S'-shapes of the nonlinear, diffraction and total force densities 
prior to and after the collision of the soliton with the right boundary; this col- 
lision is accompanied by large increases in the magnitudes of these densities. 
During the collision, the nonlinear force density shows a relative maximum 
near to the right boundary, whereas the diffraction and total force densities 
exhibit extrema. Figure 3 also indicates that, during the collision with the 
right boundary, the diffraction force exceeds the magnitude of the nonlinear 
one at the boundary, whereas the latter is larger than the former near to, but 
away from, the boundary. The momentum density shown in Figure 3 changes 
sign upon the collision of the soliton with the boundary, and its value at the 
boundary is zero. 

Comparisons amongst the momentum densities for the periodic, Dirichlet 
and Neumann boundary conditions indicate that the soliton penetrates into 
the boundary in the periodic and Neumann cases, i.e., u(L,t) ^ 0, whereas 
it docs not penetrate into the boundary for the Dirichlet case; the momen- 
tum preserves its sign in the periodic case owing to the transparency of the 
boundary and has the same shape for the Dirichlet and Neumann boundary 
conditions. 

Since the nonlinear Schrodinger equation with Robin boundary conditions 
is not a symmetric problem, the interaction of a soliton with the left boundary 
is expected to be different from that with the right one. For this reason, the 
soliton interactions with both boundaries must be considered as shown in 
Figures 4 and 5 which correspond to the first and second collisions, i.e., 
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a collision with the right boundary followed by another one with the left 
boundary. Figure 4 indicates the S'-shapes of the nonlinear, diffraction and 
total force densities prior to and after the collision of the soliton with the 
right boundary. The nonlinear force density exhibits a relative maximum and 
a relative minimum near to, but away from, the boundary where its value 
is not zero. Similar trends are observed in the diffraction force density until 
its relative maximum becomes an extremum at the boundary. Figure 4 also 
shows that the nonlinear force is greater than the diffraction one except very 
near to the boundary during the collision process. The momentum density 
presented in Figure 4 has the same shape as those for the Dirichlet and 
Neumann boundary conditions, except that its value at the right boundary 
is different from zero. 

Figure 5 indicates the lack of symmetry of the nonlinear Schrodinger 
equation subject to Robin boundary conditions. In particular, the nonlinear 
force density maintains its S-shape prior to and after the collision of the 
soliton with the left boundary; however, compared with the collision with 
the right boundary, this density exhibits a plateau before reaching a negative 
extremum at the boundary. Figure 5 also shows that the diffraction force 
preserves its S-shape prior to and after the collision of the soliton with the left 
boundary, and that it exhibits a plateau near to and a positive extremum 
at the left boundary. The total force density indicates that the nonlinear 
force is larger than the diffraction one, while the momentum changes sign 
upon the collision of the soliton with the left boundary. The isocontours 
presented in Figures 4 and 5 indicate that the soliton becomes closer to the 
right boundary than to the left one. 
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In Figures 6-10, the negative values of the kinetic, potential and total 
energy, and of the momentum are presented as functions of time. Note that 
these quantities are the integrals, from the left to the right boundary, of their 
corresponding densities, that the kinetic energy has been calculated from the 
real part of its density, and that the momentum has been obtained from the 
absolute value of its density. 

The results presented in Figure 6 indicate that the kinetic and potential 
energies are nearly constant and have opposite signs for the two-point, peri- 
odic, initial- value problem. The total energy, i.e., the sum of the potential 
and kinetic energies, is almost constant except for a minimum at about t=50 
which corresponds to the time at which the soliton crosses the boundary. 
The momentum presented in Figure 6 is also nearly constant. 

Figure 7 corresponds to Dirichlet boundary conditions and shows that 
the kinetic energy is larger than the potential one which is negative; the 
total energy is positive and exhibits a maximum at about t=50 which cor- 
responds to the collision of the soliton with the right boundary. The results 
presented in Figure 7 clearly indicate that the collision of the soliton with 
a Dirichlet boundary is accompanied by an increase (decrease) in energy as 
the soliton approaches (rebounds from) the boundary, and that the energies 
after rebound are the same as those prior to the collision. The momentum 
presented in Figure 7 indicates the deceleration (acceleration) of the soliton 
as it approaches (recedes from) the right boundary; however, these deceler- 
ation and acceleration processes do not cause large variations in the soliton 
momentum. 

The results presented in Figure 8 correspond to Neumann boundary con- 
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ditions and exhibit similar trends to those presented in Figure 7 except for 
the small, relative minima that surround the relative maxima of the kinetic, 
potential and total energies. Figure 8 also shows that the collision of the soli- 
ton with a Neumann boundary produces larger increases in energy in smaller 
intervals of time than that with a Dirichlet boundary. Large increases are also 
observed in the momentum as the soliton collides with the right boundary. 
However, both the energies and the momentum recover their values prior to 
the collision after the soliton rebounding process is completed. 

The energies presented in Figures 9 and 10 correspond to the first col- 
lision of the soliton with the right and left boundaries, respectively, for the 
Robin boundary conditions, and exhibit the same trends as, but have smaller 
magnitude than, those corresponding to the Neumann boundary conditions 
presented in Figure 8; the duration of the soliton interaction with the bound- 
aries is shorter (longer) than those of the Dirichlet (Neumann) boundary 
conditions presented in Figure 7 (8). The most noteworthy features of the 
results presented in Figures 9 and 10 are the increase in the background radi- 
ation as the soliton collides with the left boundary, the decrease (increase) in 
the momentum as the soliton collides with the right (left) boundary, and the 
relative maxima and minimum observed in the momentum upon the collision 
of the soliton with the left boundary. 

7 Conclusions 

A quantum mechanics analogy is used to determine the forces on and the 
energies of solitons governed by the nonlinear Schrodinger equation subject 
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to periodic and to homogeneous Dirichlet, Neumann and Robin boundary 
conditions in finite intervals. For all the boundary conditions considered in 
this paper, it has been shown that the the nonlinear, diffraction and total 
force densities have S-shape profiles and increase in magnitude as the soliton 
interacts with the boundaries. This interaction is characterized by an increase 
in the force densities; the largest increase corresponds to the homogeneous 
Neumann boundary conditions. 

The values of the force and energy densities after the soliton rebound 
from the boundary is completed, are same values as those prior to the inter- 
action with of the soliton with the boundary, and, except during the collision 
process, the nonlinear force density is larger than the diffraction one. During 
the collision process, the nonlinear force is larger than the diffraction one for 
the Dirichlet boundary conditions; the former is smaller than the latter at 
the boundary for the Neumann conditions. 

The lack of symmetry of the Robin boundary conditions has been illus- 
trated by the different dynamics of the soliton collisions with the right and 
left boundaries. In the collision with the right (left) boundary, the magnitude 
of the diffraction (nonlinear) density force at the boundary is larger (smaller) 
than the nonlinear (diffraction) one, and the soliton becomes closer to the 
right than to the left boundary. Since the nonlinear force is larger than the 
diffraction one when the soliton is sufficiently far away from the boundaries, 
it may be concluded that self-focusing effects are more important as regard 
the quantum mechanics analogy developed in this paper. 

The kinetic, potential and total energies of the soliton decrease as the 
softness of the boundary conditions is increased, i.e., they are largest (small- 
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est) for Neumann (Dirichlet) boundary conditions, and are nearly constant 
for the periodic boundary conditions. 
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Figure Captions 

1. Nonlinear (fn), diffraction (Re(fd)) and total (ftotal) force densities and 
momentum density (p) for periodic boundary conditions. 

2. Nonlinear (fn), diffraction (Re(fd)) and total (ftotal) force densities and 
momentum density (p) for Dirichlet boundary conditions. 

3. Nonlinear (fn), diffraction (Re(fd)) and total (ftotal) force densities and 
momentum density (p) for Neumann boundary conditions. 

4. Nonlinear (fn), diffraction (Re(fd)) and total (ftotal) force densities 
and momentum density (p) for Robin boundary conditions and first 
collision with the right boundary. 

5. Nonlinear (fn), diffraction (Re(fd)) and total (ftotal) force densities 
and momentum density (p) for Robin boundary conditions and first 
collision with the left boundary. 

6. Kinetic, potential and total energy, and momentum for periodic bound- 
ary conditions. 

7. Kinetic, potential and total energy, and momentum for Dirichlet bound- 
ary conditions. 

8. Kinetic, potential and total energy, and momentum for Neumann bound- 
ary conditions. 

9. Kinetic, potential and total energy, and momentum for Robin boundary 
conditions and first collision with the right boundary. 
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Kinetic, potential and total energy, and momentum for Robin boundary 
conditions and first collision with the left boundary. 
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Figure 1. Nonlinear (f n )> diffraction (Re{/d)) and total (/-total) force densities and 
mo merit urn density (p) For periodic boundary conditions. 
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Figure 2. Nonlinear (/„), diffraction (Rjs{J d )) and total (/^total) force densities and 
momentum density (p) for Dirichlet boundary conditions. 
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Figure 3. Nonlinear (/„), diffraction (B£{f d )) and total [/-total) fores densities and 
momentum density (p) for Neumann boundary conditions. 
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Figure 4. Nonlinear (/«), diffraction and total (/-r.ota.1) farce, densities and 

momentum density (p) for Robin boundary conditions and first collision with the 
right boundary. 
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Figure 5 Nonlinear (/u) s di fir act ion (Re(/d)) and total (/-total) force densities and 
momentum density (p) for Robin boundary conditions and first collision with the left 
boundary. 
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Figure 6- Kinetic, potential and total energy, and momentum for periodic boundary 
conditions. 
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Figure 7. Kinetic, potential and tptal energy, and momentum for Dirichlet boundary 
conditions. 
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Figure 9. Kinetic, potential and total energy, and momentum for Robin boundary 
conditions and first collision with the right boundary. 
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Figure 10. Kinetic, potential and total energy, and momentum for Robin boundary 
conditions and first collision with the left boundary. 
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